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Abstract 

The universal density-force relation is analyzed and the correspondent universal amplitude ratio 
B rea i is obtained using the massive field theory approach in fixed space dimensions d = 3 up to 
one-loop order. The layer monomer density profiles of ideal chains and real polymer chains with 
excluded volume interaction in a good solvent between two parallel repulsive walls, one repulsive 
and one inert wall are obtained. Besides, taking into account the Derjaguin approximation the 
layer monomer density profiles for dilute polymer solution confined in semi-infinite space containing 
mesoscopic spherical particle of big radius are calculated. The last mentioned situation is analyzed 
for both cases when wall and particle are repulsive and for the mixed case of repulsive wall and 
inert particle. The obtained results are in good agreement with previous theoretical results and 
with the results of Monte Carlo simulations. 
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The investigation of polymer solution near the surfaces and in film geometries or meso- 
scopic particles dissolved in the solution is a task of great interest from theoretical and 
application point of view. 

The monomer density profiles of dilute polymer solution bounded by a planar repulsive 
wall has a depletion region of mesoscopic width of order of the coil size R x ~ N u (where N 
is the number of monomers per chain), and for the distances z from the wall that are small 
compared to this width but bigger than microscopic lengths of monomer size / the profile 
increases as 

p(z) ~ z~ 

with Flory exponent v {y is 1/2 for ideal polymer chains and v ~ 0.588 for polymer chains 
with excluded volume interaction (EVI)). This remarkable theoretical predictions was pro- 
posed by Joanny, Leibler and de Gennes [1]. They also mentioned that the monomer density 
close to the wall is proportional to the force per unit area which the polymer solution exerts 
on the wall. But, for the first time a complete quantitative expression for the universal 
density-force relation was obtained by Eisenriegler on the basis of e- expansion up to first 
order in [2]. As was mentioned in [2], the correspondent density-force relations with the same 
universal amplitude B are valid for the different cases: 1) a single polymer chain with one 
end (or both ends) fixed in the half space bounded by the wall; 2) a single chain trapped in 
the slit of two walls; 3) for the case of dilute and semi-dilute solution of free polymer chains 
in a half space; 4) for the case of polymer chain in a half space containing a mesoscopic 
particle of arbitrary shape. The verification of the universal density-force relation was per- 
formed by simulation techniques using an off-lattice bead-spring model of a polymer chain 
trapped between two parallel repulsive walls [3] and by the lattice Monte Carlo algorithm 
on a regular cubic lattice in three dimensions [4]. Unfortunately, the obtained results of 
Monte Carlo simulations in [3] are much higher B MC > 2.48 than theoretical predictions 
B e 1.85 obtained in [2]. As it was mentioned in [3], there is a systematic decrease of B e ff 
with increasing the distance L between the walls. In accordance with it some rough linear 
extrapolation with plotting B e ff versus l/y/L which yields an extrapolated result B ~ 1.4 
was performed in [3]. Recent numerical results based on the lattice Monte Carlo algorithm 
on a regular cubic lattice gives the value B ps 1.70 ± 0.08 which is smaller than previous 
theoretical predictions. This all indicates that the above mentioned task still have a lot of 
open questions and the present paper tries to give answers for some of them. 
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The present paper is devoted to investigation of the universal density-force relation and 
calculation of the universal amplitude B in the framework of the massive field theory ap- 
proach. The massive field theory approach gives better agreement with the experimental 
data and the results of Monte Carlo calculations as it was shown in the case of infinite [5, 6], 
semi-infinite [7] systems, and specially in the case of dilute polymer solutions in semi-infinite 
geometry [8] and confined geometry [9]. As was mentioned above, the knowledge of the uni- 
versal amplitude B allows to obtain the monomer density profiles for whole class of different 
systems (cases l)-4)). Besides, it should be mentioned that the universal density- force rela- 
tion is valid not only for the case of polymer chains trapped in between two repulsive walls 
but also for the mixed case of one repulsive and one inert wall. 

We consider a dilute polymer solution, where different polymer chains do not overlap and 
the behavior of such polymer solution can be described by a single polymer chain. As it 
is known, the single polymer chain can be modeled by the model of random walk and this 
corresponds to the ideal chain at ^-solvent or self-avoiding walk for the real polymer chain 
with EVI for temperatures above the #-point. Taking into account the polymer-magnet 
analogy developed by [10], their scaling properties in the limit of an infinite number of steps 
N may be derived by a formal n — > limit of the field theoretical 4 0(n)- vector model 
at its critical point. It should be mentioned that 1/N plays the role of a critical parameter 
analogous to the reduced critical temperature in magnetic systems. The role of a second 
critical parameter plays the deviation from the adsorption threshold (c oc (T—T a )/T a ) (where 
T a is adsorption temperature). The value c corresponds to the adsorption energy divided by 
ksT (or the surface enhancement constant in field theoretical treatment). The adsorption 
threshold for long-flexible polymer chains, where 1/N — » and c — > is a multicritical 
phenomenon. 

In order to obtain the universal amplitude B in universal monomer-density relation let's 
consider the single polymer chain with one end fixed fluctuating near the repulsive wall such 
that (I « z « z « R x ). 

The correspondent layer monomer densities p\{z) defined by [2] is: 

(R ) 1/u 

Px (z)dz= [ -f^ n dN x (z), (1) 

where dN\(z) is the number of monomers in the layer between z and z + dz, R x is the 
projection of the end to end distance R onto the direction of x axis. Besides, p\ is obtained 
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from monomer density p(r, z) after integration over the d — 1 components parallel to the 
wall. The scaling dimension of p(r, z) is l 1 / u ~ d and equals the ordinary dimensions of the 
quantity 

*(*) = l|^ 2 (*)> (2) 
where L is the conjugate Laplace variable which has the dimension of length squared 
and is proportional to the total number of monomers N of the polymer chain. Following 
the description of the problem as given in [2], the monomer density in this case is 

XC < 0(x)0(x') > w 

in the limit n — > 0. The average < • > w in (3) denotes a statistical average for a Ginzburg- 
Landau field in semi-infinite geometry. The dot in (3) means a cumulant average. The 
correspondent Ginzburg-Landau Hamiltonian describing the system in semi-infinite (j — 1) 
or confined geometry of two walls (j = 1, 2) is: 

+ E % / rfd ~V 2 , (4) 

— * 

where 0(x) is an n- vector field with the components <f>i(x), i = 1, n and x = (r, z), fj, 
is the "bare mass", v is the bare coupling constant which characterizes the strength of the 
excluded volume interaction (EVI). The surfaces of the system is characterized by a certain 
surface enhancement constant Cj , where j = 1,2. The interaction between the polymer 
chain and the walls is implemented by the different boundary conditions. 

In the case of two repulsive walls (where the segment partition function and thus the 
partition function for the whole polymer chain tends to as any segment approaches the 
surface of the walls) Dirichlet-Dirichlet boundary conditions (D-D b.c.) takes place: 

C\ — > +oo, c 2 — > +oo or 0(r, 0) = 0(r, L) = 0, (5) 

and for the mixed case of one repulsive and one inert wall Dirichlet-Neumann boundary 
conditions(D-N b.c.) are: 

Cl ^+cx), c 2 = or 0(r,O) = O, d ^ Z \ z=L = 0. (6) 



Taking into account that J d d x < p(x) >= R l J u the property 

J d d HC < *(x) • 0(x)0(x') > w = R l J v XC < 0(x)0(x') > w , (7) 

takes place. It should be mentioned that near the repulsive wall the short-distance expansion 
of <fr 2 takes place and it assumes [11-14] 

tt(M) -+ Bz 1 '" ^- 1 ^ 2 , (8) 

for the distances I « z, where I is monomer size. The surface operator with 
&± — 2 " > \s=o is the component of the stress tensor perpendicular to the walls. Taking 
into account the shift identity [2, 15, 16] for the case of semi-infinite geometry 

J d d -H < 2^ ■ 0(x)0(x') >»= (§~ z + ^) < 0(x)0(x') > w (9) 

and integrating it over J d d x. f for the layer monomer density Eq.(l) in accordance with 
Eqs.(3), (8) the universal density-force relation can be obtained 

< p x (~z) >= Bz^-J-, (10) 

where 

^ = A l n [J£ J rfV < 0(x)0(x') >J (11) 

is the force per area that the polymer chains exert on the wall. It should be mentioned 
that the density-force relation (10) takes place for the distances / << z « R x , and B 
is universal amplitude. Following the scheme of obtaining the universal amplitude B as it 
was proposed in [13, 16] the correspondent calculations based on the massive field theory 
approach in fixed d = 3 dimensions were performed. Thus, for the layer monomer density 
of ideal polymer chain takes place: 



< Rl >V*» XCG{- z, z) / °° dz'G(; ~z, z>) 



z2 



Z 



< > 0= Icfr z G ( ; ^) = 2 ~ (12 > 

Taking into account the value of f/ksT ps \ jz for the force exerted by ideal polymer chain 
on the wall from Eq.(10) the universal amplitude B can be obtained: 

B id = 2. (13) 
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The case of real polymer chain is more complicated, because EVI with nonequal to zero 
the bare coupling constant t> should be taken into account. Taking into account Eq.(lO) 
and Eqs.(3),(ll) after renormalization of the mass 

nl = H 2 -jJM + 0(v>), (14) 

where Ji(/Jo) = / d 2K 9 with K q = q 2 + /Xq, the renormalization of the coupling 

constant v = vji and including the correspondent UV-finite renormalization factors (see 
[17]) in the limit n — > up to one- loop order for the real polymer chain with EVI the 
universal amplitude B rea i can be obtained : 

£ = £^ = 2(1-^(1 + !^ - 7 *)), (is) 

where ^ E = 0.577 is Euler's constant. Here we took into account that v = |(1 + |) + 
0(v 2 ) and the following definition v = b n (d)v was introduced with b n (d) = ^§ ^r[l/2] • ^ ne 
correspondent fixed point is equal: v* — 1. At d = 3 Eq.(15) leads to: 

Srea/ ~ 1-62. (16) 

The obtained result Eq.(16) is in agreement with the result obtained by Eisenriegler [2] for 
real polymer chains in the framework of e-expansion at d — 3: B — B^ al = 2(1 — be) m 1.85, 
where e = 4 — d and b = (In 2 + ^ E — 2/3)/8. As it is easy to see, the result obtained in the 
framework of the massive field theory approach is slightly smaller than result of e-expansion 
[2, 13] and is in agreement with the numerical result B rea i 1.70 ± 0.08 of Monte Carlo 
simulations [4] and with the result of rough linea extrapolation B extr ~ 1.4 obtained on the 
basis of plotting B e ff versus 1/y/Z in [3]. 

As it was shown in [2], the density-force relation (10) is also valid for the case of dilute and 
monodisperse solution of free chains in semi-infinite space. The polymer density far from 
the wall is fixed and the pressure on the wall f/A (where A is the area of the wall) equals 
the pressure in the bulk risksT . Thus, the monomer density of dilute polymer solution of 
free chains in semi-infinite space according to Eq.(10) is 

<p(x)> f =Bz 1 '"n B , (17) 

where n# is the polymer density in bulk far from the wall. 

In the case of a spherical particle with radius R much larger than the distance of its 
closest point a to the surface and much larger than radius of gyration R g the Derjaguin 
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approximation, which describes the sphere by a superposition of fringes with local distance 
from the wall £( r ||) = a + r |/(2R-) ? should be applied [18]. 

Taking into account the depletion interaction potential between the particle and the wall 
which we obtained in [9] (see Eq.(7.12)): 



oo 

n B k B T 



J dye(y), (18) 



a/R x 

and the correspondent scaling function for the free energy of interaction for the slit geometry 
Q(y) (see Eqs.(5.6),(5.8) and Eqs.(6.4),(6.8) in [9]) the layer monomer density of dilute 
polymer solution in semi-infinite space containing spherical particle of big radius for A = 1 
and ub = I /{LA) can be obtained: 

< Px(z) > WP = ^(1 - 2nRa 2 Q(^)). (19) 

Li lt x 

As it is easy to see from Eq.(19), the layer monomer density depends not only on R x , but 
also on the shape of the mesoscopic particle and its distance from the wall. 

The density-force relation in the case of single polymer chain trapped in the slit of two 
walls can be shown in the same way, because the Eq.(7) and Eq.(8) take place not only for the 
averages of the type <•>„,, but also for the averages < • > ww . More detailed calculations 
can be found in [17]. As it was mentioned in [2], the monomer density in the case of dilute 
polymer solution between two repulsive walls has a maximum in the center z = L/2 of 
the slit. If the distant wall at z = L is inert or in another words is at the adsorption 
threshold, the density-force relation Eq.(lO) again takes place with the same values of B 
(see Eqs.(13),(16)), as it was mentioned by Eisenrieglier [2]. In this last mentioned case the 
polymer chain prefers the distant inert wall and the monomer density maximum is at z = L. 
The results of calculations of the layer monomer density profiles for the case of ideal polymer 
chain and real polymer chain with EVI confined between two repulsive walls (D-D b.c), one 
repulsive and one inert wall (D-N b.c.) in accordance with Eq.(lO) and taking into account 
the correspondent values of B^ai and B rea \ (see Eqs.(13), (16)) are presented on Fig.l and 
Fig. 2, respectively. Besides, Fig.l and Fig. 2 present results for the case of ideal and real 
polymer chains in semi-infinite space containing spherical particle of big radius. The last 
mentioned situation is analyzed for both cases when wall and particle are repulsive and for 
the mixed case of repulsive wall and inert particle. 
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The obtained results (see Fig.l) indicate that the layer monomer density profiles for 
ideal polymer chains are weaker then for real polymer chains with EVI in the case of two 
repulsive walls (or D-D b.c). Completely opposite behavior of monomer density profiles is 
observed for the case of one repulsive and one inert wall (or D-N b.c), as it is easy to see 
from Fig. 2. Besides, the layer monomer densities for curvature surfaces are smaller then for 
planar surfaces. 

In order to test the reliability of the obtained analytical results it will be interesting to 
compare them with the recent results of Monte Carlo calculations obtained by [4] for the 
single polymer chain trapped inside the slit of two repulsive walls. 

In Ref. [4] the lattice Monte Carlo algorithm on a regular cubic lattice in d = 3 dimen- 
sions, with D lattice units in z-direction and impenetrable boundaries was applied (L = uD 
with u denoting the lattice spacing). The other directions obeyed periodic boundary con- 
ditions. The correspondent reduced force in accordance with ([4]) can be written in the 
form: 

2 VF(#V\ (20) 



k B T vfiooL \R X 
where parameter a is a universal amplitude, ^ is the critical fugacity per monomer. Taking 
into account Eq.(10) and the value B^eai from (13), the monomer density near the wall for 
ordinary RW is scaled as: 

2ttV ( L\- 2 

< px(z) > G = —r- bH » ( 21 ) 



L \R r 

where the universal amplitude a for the case of ideal chains was found as a 0.2741(2) 
(see Ref. [4] ), which is very close to the exact value, computed analytically in [2] and equal 
to a = For the ideal chain takes place: v = 0.5, Xd — 1 and Poo = f ■ 111 Fig.3 this 
asymptotic behaviour for narrow slits is clearly recovered by our results for ideal chains, 
where the narrow slit limit is valid. 

In the case of SAW in Ref. [4] the value a « 0.448 ± 0.005 was obtained. Taking into 
account the values for v « 0.588, x\ ~ 0.960 (e.g. [19]) and ^ w 0.2135, B real = 1.70 ±0.08 
(see Ref. [4]), the correspondent monomer density for SAW can be written as 

15.44 z 1 - 7 ( L V 1 ' 7 
< Px(z) >evi~ z ( 7T ) • ( 22 ) 



R 



The result Eq.(22) is presented in Fig.3 and compared to our theoretical results for a trapped 
chain with EVI, which are valid for the wide slit regime y > 1. In order to compare the 
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obtained analytical results with recent results of MC simulation we performed extrapolation 
of the MC results for a and fi^, which are valid for the narrow slit regime to the region of 
y > 1. As it possible to see from Fig.3, the results Eq.(21) and Eq.(22) correspond very 
well to our theoretical predictions in the wide slit limit. One of the possible reasons for 
remaining deviations with the results of Ref. [4] is that the chain in the MC simulation is 
too short in order to compare with the results of the field-theoretical RG group approach. 
Unfortunately, at the moment no simulations concerning one repulsive and one inert wall 
exist. 
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FIG. 1: The dimensionless value of the layer monomer density < p\{z) > L 1 ~ 1 / u profiles for: ideal 
and real polymer chains in the case of two repulsive walls (D-D b.c)(ww) with y = 10, and dilute 
polymer solution (R x = 0.1pm) of ideal and real polymer chains in semi-infinite space containing 
repulsive spherical particle of big radius R = 10.0pm on the distance a = 1.0pm from the repulsive 
wall (a = L) (pw). The maximum is at L/2. 
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FIG. 2: The dimensionless value of the layer monomer density < p\(z) > L 1 ~ 1 / u profiles for: ideal 
and real polymer chains in the case of one repulsive one inert wall (D-N b.c.) (ww) with y = 10, 
and dilute polymer solution (R x = 0.1pm) of ideal and real polymer chains in semi-infinite space 
containing inert spherical particle of big radius R = 10.0pm on the distance a = 1.0 pm from the 
repulsive wall (a = L) (pw). The maximum is near the distant inert wall. 




FIG. 3: Comparison of theoretical results with Monte Carlo simulations for the layer monomer 
density < p\{z) > profiles for: ideal chain and real polymer chain with EVI (MFT) in the case 
of two repulsive walls (D-D b.c). RW(MC) and SAW(MC) are the results obtained in Ref.[4] 
for random walks and self avoiding walks. The comparison of results are performed for the case: 
D = 80, N = 3500 and y rj 2.34 (RW); D = 120, N = 3500 and y ps 1.73 (SAW) . 
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